Effect of Ordering on Spinodal Decomposition of Liquid-Crystal/Polymer Mixtures 
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Partially phase-separated liquid-crystal/polymer dispersions display highly fibrillar domain mor- 
phologies that are dramatically different from the typical structures found in isotropic mixtures. 
To explain this, we numerically explore the coupling between phase ordering and phase separation 
kinetics in model two-dimensional fluid mixtures phase separating into a nematic phase, rich in 
liquid crystal, coexisting with an isotropic phase, rich in polymer. We find that phase ordering can 
lead to fibrillar networks of the minority polymer-rich phase. 
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Mixtures of liquid crystals with a small amount of poly- 
mer (polymer-stabilized liquid crystals or PSLC's) 
show promise for electro-optic devices such as light shut- 
ters and displays Q-Q], because the polymer tends to 
form a network that aligns the liquid crystal Since 
polymers and liquid crystals tend to be immiscible, the 
dispersions are prepared by mixing a small amount of 
miscible monomer with the liquid crystal and photopoly- 
merizing. As the polymers grow, the system phase sep- 
arates into an ordered phase rich in liquid crystal and 
an isotropic phase rich in polymer. Long before the 
system reaches equilibrium, however, the polymeriza- 
tion "freezes" the mixture into a crosslinked network of 
polymer-rich domains. Thus, the fabrication of PSLC's 
involves interplay among three kinetic processes: poly- 
merization, phase separation, and phase ordering. De- 
pending on the time scales that control these processes, 
a rich variety of morphologies have been observed [p|-^2| . 
Because of the number of nonequilibrium processes in- 
volved, however, there is little theoretical understanding 
of the factors that control the domain morphology. In 
this Letter, we focus on the interplay between phase sep- 
aration (PS) and phase ordering (PO) kinetics in mix- 
tures of short, rigid polymers (rods) and long, flexible 
polymers (coils), as a first step towards rational design 
and control of the network morphology. 

It is well known that thermodynamic factors such 
as the anisotropy of the isotropic/nematic interfacial 
tension can influence domain morphology, leading to 
anisotropic domain shapes. However, there are also ki- 
netic factors that control domain morphology, such as 
the anisotropic diffusion coefficient of a rod. To capture 
these thermodynamic and kinetic effects, we use a Cahn- 
Hilliard framework that allows composition and orienta- 
tional density to evolve in a coupled fashion as functions 
of position and time following a temperature quench [fl3| . 
In contrast to earlier studies that treat orientational den- 



sity as a scalar order parameter p^Jl5[ |, this framework 
includes the orientational density's second-order tensorial 
nature [^6|. Although it is instructive to study the case 
of two coupled scalar order parameters (Model C pTf), 
a scalar cannot capture the direction of nematic order. 
Because a vector does not have head/tail symmetry, it 
is crucial to retain the tensor order parameter to obtain 
domain anisotropy |l8| ]. 

To assess the effects of phase ordering, we study two 
systems. The first (denoted RC) is a mixture of rods 
and coils. The second (denoted CC) is a polymer blend, 
identical to RC except that the rods are replaced by flex- 
ible chains of the same length that do not align. In 
both cases, we solve the linearized coupled partial dif- 
ferential equations of motion analytically, and the non- 
linear equations numerically into the late-time regime, 
after quenching the isotropic, homogeneous mixture into 
the coexistence region. We find that phase ordering dra- 
matically affects morphology, giving rise, for example, to 
interconnected networks or elongated domains depend- 
ing on whether phase ordering or phase separation is the 
initially dominant process. 

We study two-dimensional, incompressible mixtures of 
chain molecules made up of monomers of the same size. 
In RC, we have short rods that are iV^ = 10 monomers 
long, and coils that are Nb — 100 monomers long. The 
phase behavior of the system is governed by a bulk free 
energy |l^,^l -Frc that couples the area fraction of rods, 
(f>, to the orientational density of the rods, S. This free 
energy consists of two parts: the Flory-Huggins free en- 
ergy, which governs phase separation, and the Landau-de 
Gennes free energy, which governs the isotropic/nematic 
transition. The free energy was calculated within the 
random phase approximation |2(J]. It depends on the 
rod and coil lengths, and contains only two free param- 
eters. These are the two interaction parameters, namely 
the Flory parameter, \, and the Maier-Saupe parameter, 
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w, which control the isotropic and anisotropic monomer- 
monomer interactions, respectively. Both of these param- 
eters depend on temperature with cntropic and enthalpic 
contributions: \ — Xo + Xi/T and w = wq + Wi/T. We 
have chosen the parameters \% an d w i so that the critical 
temperature for phase separation according to the Flory- 
Huggins free energy, T c , and the isotropic/nematic tran- 
sition temperature of the pure rod system, T n i , are both 
of order unity in our (arbitrary) units, with T n i higher 
than T c . We use xo = 0.055, xi = 0.036, w = 0.55, and 
wt = 0.26, which yield T c = 1.14 and T ni = 1.3. The 
resulting RC phase diagram is shown in Fig. |l|. 
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FIG. 1. Phase diagram of two-dimensional rod/coil blend 
for the parameters discussed in the text. T is the temperature 
in arbitrary units chosen so that the IN transition tempera- 
ture is of order unity, and 4> is the local area fraction of rods. 
The heavy solid line denotes the coexistence curve. At low 
temperatures, an isotropic phase rich in coils coexists with a 
nematic phase rich in rods. 

Below a triple point, the system demixes into an 
isotropic phase rich in coils and a nematic phase rich 
in rods. Above the triple point, there is an isotropic- 
isotropic coexistence region and an isotropic-nematic 
(IN) coexistence region. The IN coexistence region ends 
in a critical point that is connected to the IN transition 
of a pure rod system by a second-order transition line. 
Note that in three dimensions, the IN transition is first 
order, so the phase diagram would be different. Below 
the triple point, however, we find that the IN coexistence 
curves are qualitatively identical in two and three dimen- 
sions. The dashed line represents the spinodal for phase 
separation (PS), while the dotted line is the spinodal for 
phase ordering (PO). 

In the CC system, we also use — 10 and Ng = 100 
monomers. Because the "rods" are now flexible, the free 
energy Fee is given by the Flory-Huggins free energy 
only, without the Maier-Saupe contribution [|T). The 
parameter x is exactly the same as in the RC system. 
Thus, CC is identical to RC except that there is no phase 
ordering. 



We study the morphology following four quenches from 
the isotropic, homogeneous phase into the IN coexistence 
region below the triple point, to the four points marked 
A-D in the RC phase diagram. We use a simplified ver- 
sion of the equations of motion derived within the dy- 
namical random phase approximation |l3| 

at = r ^ v { ) 

dt ~ 33 5S kl { ) 

where and are calculated Onsager coefficients 
that depend on the single chain and single rod diffusion 
coefficients; expressions for these coefficients are provided 
in Ref. [|l3|. The free energy functional -F(</>, S) consists 
of the bulk free energy and a nonlocal free energy that 
controls the cost of gradients in composition and orien- 
tational density |20|] . Note that the cost of gradients in 
orientational density are determined by the Frank elastic 
constants that characterize nematic elasticity. For CC, 
the time evolution is given by Eq. [I] alone, since S = 
and SFcc/SS^ = 0. We emphasize that for both sys- 
tems, the equations of motion are not phenomenological; 
they are derived from microscopic models. We solve the 
discretized equations on a square lattice of size 150 2 using 
a variable timestep Runge-Kutta numerical integration 
scheme. We have also varied the number of lattice points 
and the size of the mesh to verify that our observations 
do not depend on system size or discretization. 

Quench A: Here we quench to a point where RC is ini- 
tially unstable with respect to PS, but metastable with 
respect to PO. Thus, the point marked A in Fig. [I] lies 
below the spinodal for PS and above the one for PO. For 
CC, where phase ordering is absent, this quench leads 
to circular droplets rich in the shorter coils, in a matrix 
rich in the longer coils. RC initially displays almost iden- 
tical behavior with circular droplets rich in rods. Once 
the concentration of rods in the droplets is comparable 
to the concentration at the IN spinodal, however, the 
rods in the droplets begin to order. At this time, the 
droplets abruptly expel more coils because the coils are 
less soluble when the rods are ordered. In addition, the 
droplets must develop defects because the rods want to 
be parallel to each other and parallel to the droplet in- 
terface. As the rods order, a pair of defects forms inside 
each droplet and separates, with the two defects mov- 
ing along the director in opposite directions towards the 
edge of the droplet. The magnitude of orientational or- 
der is lower at the edges than in the center; thus, the 
defects migrate towards the edges to lower their energy. 
Since the overall system is isotropic, the long axis of each 
droplet is randomly oriented. Note that it is essential to 
retain the tensorial nature of the orientational order pa- 
rameter in order to obtain randomly-oriented elongated 
domains. Fig. p] shows a snapshot of the RC system at 
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a time at which there is significant phase separation and 
ordering. 




FIG. 2. Late-time snapshot of the local composition for the 
RC system when quenched to A, as indicated in Fig. |l| We 
use a gray scale where black corresponds to pure coils and 
white corresponds to pure rods. Arrows represent the local 
nematic director (shown for two droplets only). 

Quench B: For this quench, CC is metastable with re- 
spect to PS, so isolated droplets rich in the long coils 
would form by nucleation if we had included thermal 
noise. Since we have not included thermal noise, the CC 
system does not evolve. In this region RC is also initially 
metastable with respect to PS, but it is unstable with 
respect to PO. The instability towards orientational or- 
dering eventually drives the system to phase separate be- 
cause the two order parameters are coupled. This makes 
physical sense: once the rods are strongly aligned, they 
expel the coils into isolated droplets. Although there is 
no orientational order within the coil-rich droplets, the 
droplets are anisotropic because of the nematic elasticity 
in the surrounding rod-rich matrix. 

The appearance of droplets even when the system is 
only metastable to phase separation is reminiscent of 
recent theoretical work on protein and polymer crys- 
tallization, which shows that a quench below the spin- 
odal for one order parameter (density) can promote do- 
main growth in another order parameter (crystallization) 
p2|p3j . The advantage of our system is that we can cal- 
culate the coupled equations of motion for our two order 
parameters ^from a microscopic model of rods and coils. 

Quench C: In this quench, RC is initially unstable both 
to PS and to PO, but is more unstable with respect to 
PS. In both RC and CC, the system forms a bicontinuous 
network that rapidly breaks up into droplets rich in long 
coils. In RC, however, defects in the surrounding rod- 



rich domains give rise to droplets that are noncircular, 
as shown in Fig. |3[ Furthermore, the coupling of PS and 
PO leads to a faster onset of phase separation in RC than 
in CC, because the rod-rich regions tend to expel coils as 
the rods order. 




FIG. 3. Snapshot of the local composition for the RC sys- 
tem at a late time following a quench to C. Black corresponds 
to pure coils and white corresponds to pure rods. Arrows il- 
lustrate the local nematic director, and a defect is marked 
with a circle. 

Quench D: As in Quench C, RC is initially unstable 
both to PS and to PO. However, in this quench it is more 
unstable with respect to PO. Once the degree of order is 
significant, then phase separation begins. Because of PO, 
PS begins much earlier in RC than in CC; by the time CC 
displays any noticeable phase separation, the concentra- 
tion difference between the two phases in RC has already 
reached its equilibrium value. Fig. || presents snapshots 
of the systems as they evolve in time. In CC, an intercon- 
nected network rich in the longer coils initially appears. 
However, it quickly breaks up into droplets that become 
increasingly circular over time. In RC, on the other hand, 
small highly-ordered rod-rich droplets initially form in a 
coil-rich matrix. As these rod-rich droplets grow, the sur- 
rounding coil-rich region shrinks into a network of inter- 
connected domains that are strikingly fibrillar, despite 
the fact that the coil-rich phase is the minority phase. 
This is in contrast to systems with only a compositional 
order parameter (Model B fT^] ), where only the major- 
ity phase can form networks. Eventually, this network 
breaks up to form coil-rich droplets; again, these are non- 
circular because of defects in the surrounding rod-rich 
regions. However, in experiments the crosslinking of the 
polymer network may arrest the phase separation at the 
fibrillar network stage. 
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black corresponds to pure long coils in CC and to pure coils 
in RC, while white corresponds to pure short coils in CC 
and pure rods in RC. Adjacent snapshots for the two sys- 
tems have comparable degrees of phase separation. Note that 
phase ordering in RC drives the system to a fibrillar network 
morphology of coil-rich domains at intermediate times. 

In summary, we find that phase ordering can signif- 
icantly influence domain morphology, especially when 
(1) the system is unstable to both phase separation and 
phase ordering and (2) orientational order is significant 
at the onset of phase separation. This is consistent with 
morphologies observed in PSLC's, where fibrillar net- 
works of the minority polymer-rich domains form when 
the system is initially ordered [p|-pd|. Our results show 
that we can learn much about domain morphology from 
the phase diagram and spinodal lines, and that the cou- 
pling of phase separation and ordering kinetics leads to 
rich phenomena that warrant further study. 
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